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ABSTRACT 



We study on-shell and off-shell properties of the supersymmetric sinh-Gordon 
and perturbed SUSY Yang-Lee models using the thermodynamic Bethe ansatz 
and form factors. Identifying the supersymmetric models with the Eight Ver- 
tex Free Fermion Model, we derive inversion relation for inhomogeneous transfer 
matrix and TBA equations and get correct UV results. We obtain two-point 
form factors of the trace of energy-momentum tensor using the Watson equa- 
tions and their SUSY transformations. As an application, we compute the UV 
central charge using these form factors and spectral representation of the C- 
theorem. 



1. Introduction 



For 2D integrable field theories ^'-matrices are purely elastic, all incoming 
momenta are conserved and multi-particle scattering amplitudes are factorized into 
a product of two-particle S'-matrices. These S'-matrices, in turn, should satisfy 
Yang-Baxter equations which often determine the S'-matrices completely along 
with unitarity and crossing symmetry [1] . The .S-matrix provides essential tools to 
understand 2D field theories. First of all, the ^'-matrix gives information on the UV 
behaviour of the theory by relating the Casimir energy on the cylinder to the central 
charge of the corresponding UV conformal field theory (CFT) [2] . This program 
known as thermodynamic Bethe ansatz (TBA) [3] has provided consisitency checks 
for many factorizable scattering theories either with local lagrangians or without 
them such as perturbed CFTs [4]. 

(S'-matrix plays an important role in off-shell physics as well. It can be used 
to determine off-shell quantities such as correlation functions by computing the 
matrix elements of an operator on the basis of the on-shell particles. These objects 
known as form factors (FFs) may be computed exactly using only the S'-matrices 
and particle spectrum (bound states) as input [5,6]. With exact FFs correlation 
functions are given by an infinite sum over intermediate on-shell states. This form 
factor approach has an advantage for the computation of correlation functions of 
massive integrable models that the infinite sum over all intermediate states con- 
verges very fast. For many cases, upto two-point FFs give quite accurate results on 
off-shell quantities [7-10]. Furthermore, the two-point FFs can be related to some 
exact non-perturbative informations of the underlying theories, such as the wave 
function renormalization [5,12] and the UV central charges through the spectral 
representation of the C-theorem [8,13,9]. In this sense, without complete solutions 
of the FFs one can still extract non-perturbative off-shell informations from the 
FFs. 

While the TBA analysis or the FF computation can be relatively simple for 
diagonal scattering theory, which has no mass degeneracy, non-diagonal scattering 
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theories entail much more comphcacy. By nondiagonal we mean theories with dif- 
ferent types of particles of the same mass for which the scattering of two particles 
can occur in more than one channels. Most of interesting 2D integrable field theo- 
ries such as the soliton scattering theories, theories with internal gauge symmetries, 
and supersymmetric theories belong to this class. 

For the non-diagonal theories, the equations for the TBA and FFs are expressed 
in terms of monodromy and transfer matrices. To solve the equations, one needs to 
diagonalize these matrices. It is remarkable that with some technical diffferences 
the same problem is often met in the study of the lattice models [14]. In the lattice 
model the Yang-Baxter equations are to be satisfied to construct infinite number 
of conserved charges through the commuting transfer matrices. Partition functions 
and free energies are expressed in terms of the eigenvalues of the transfer matrices. 
Due to this common feature, it is often quite useful to connect 2D field theories 
with lattice models. 

There are two types of the models in the lattice and continuum which are 
connected with each other. The first one is so-called the vertex type; the states 
are assigned on the lines which form a lattice. For the square lattice, each vertex 
consists of four lines and assigned a Boltzmann weight depending on the four 
states of the hues [14] . These hues correspond to the world hues of incoming and 
outgoing particles in the scattering theories. While some of these vertex models are 
associated with field theories with local lagrangians, there remain many vertex-type 
lattice models still to be related to 2D integrable field theories. 

The second type is the interacting- round- face models [15]. The Boltzmann 
weights are assigned on each vertex on the square lattice, depending on the heights 
of four faces. As a special case, if the heights are restricted, one obtains restricted 
solid-on-solid (RSOS) type of models. These wide class of the lattice models have 
been related to 2D CFTs. Due to the conformal invariance, the corresponding lat- 
tice models are at the criticality. Many exact results including correlation functions 
have been obtained using the CFT techniques. These identification can be con- 
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tinued in the off-critical region. Without the conformal symmetry, the off-critical 
RSOS models are associated with CFTs perturbed by relevent operators [16-19]. 
Again, .S-matrices of the perturbed CFTs are given by the Boltzmann weights of 
the RSOS models. 

The best known example is the relation between the six vertex model and the 
sine-Gordon (SG) model. The SG model has soliton and antisoliton spectrum and 
the (S"- matrix can be associated with the i?- matrix of the s\q{2), affine quantum 
group [18]. The Boltzmann weights of the six vertex model are the same as the 5"- 
matrix elements after identifying the up and down arrows assigned on each vertex 
line with the soliton and antisoliton. In addition, quantum group reduction of the 
SG model corresponds to the RSOS lattice model obtained from the six vertex 
model. The TEA analysis of these models have been done by diagonalizing the 
inhomogeneous transfer matrices of the six vertex [21,22,23] and RSOS models [24]. 

The complete FFs of the SG model have been obtained by F. Smirnov us- 
ing quantum inverse scattering methods, providing only known example with the 
complete FFs for nondiagonal theories. Based on this information, Smirnov found 
axioms for the FFs to satisfy [6]. Therefore, the problem to find complete FFs is 
reduced to solve these axioms for a given theory. However, solving these axiomatic 
equations completely is very difficult even for diagoanal scattering theories except a 
few simplest ones such as Ising, Yang-Lee, and sinh-Gordon models [8,7,9,11]. The 
problem becomes much more complicated for the nondiagonal cases. As an initial 
step to the problem, we will concentrate on two-point FFs. Two-point FFs can 
be determined relatively easily by diagonalizing ^'-matrix and evaluating the FFs 
using the Watson equations [5]. For the supersymmetric theories, details can be 
further simplified due to the SUSY relations between the FFs. As stressed before, 
the two-point FFs have many useful informations on the underlying theories. 

In this paper, we want to apply these frameworks to the N — 1 supersymmetric 
(SUSY) theories. The -S-matrices of many SUSY models have been otained. These 
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^'-matrices have the following factorized form [19,26]: 

s{e) = Ss{e)®So{e), (i.i) 

where the first factor carries the SUSY indices and commutes with the SUSY 
charges while the second one 5*0 is the S'-matrices of the models without the SUSY. 
So far, several SUSY integrable field theories and perturbed super CFTs are solved 
and their ^'-matrices are derived. Interesting aspect of the SUSY models is that 
these (S'-matrices commuting with SUSY charges are identified with Boltzmann 
weights of some lattice models. 

For example, for the = 1 SUSY CFTs perturbed by the least relevent 
operator, Sq, which commutes with SUSY charges with central extension due to the 
topological charges, is related to the RSOS weights corresponding to the tricritical 
Ising model [25]. For the N — 2 SUSY models, the first factor is identified with 
the Boltzmann weights of the six vertex model [20,22]. These relations with lattice 
models are important not only for the lattice- field theory correspondence but for 
actual solutions of the models. 

A'" = 1 SUSY sinc-Gordon (SSG) model has been solved in an unconventional 
way. Its soliton S'-matrix has been derived from the results on the perturbed super 
CFTs by the least relevent operator [19]. The SUSY part of the SSG soliton ,5- 
matrix is given by the RSOS tricritical Ising model ^'-matrix while Sq is ordinary 
sine-Gordon ^'-matrix. The -S'-matrices of the SSG bound states (breathers) have 
been derived from multi-soliton scattering amplitudes [27]. In particular, since 
the lightest bound states are forming a supermultiplet of the fundamental fields 
appearing in the SSG lagrangian, the lightest breather S'-matrix of the SSG model 
can be analytically continued to get the S'-matrix of supersymmetric sinh-Gordon 
(SShG) model. This S'-matrix is identical to the one derived first by Shankar and 
Witten by explicitly requiring the commutativity with SUSY charges [28] . Besides, 
the SSG model with only the lightest breather in the spectrum can be understood 
as perturbed super CFTs, the SUSY Yang- Lee (SYL) model [26,27]; the simplest 
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nonunitary super CFT perturbed by the least relevent operator. This model in- 
cludes only one supermultiplet of on-shell states and the S'-matrix is identical with 
that of the SShG model. This ^-matrix is our starting point. 

These models with N — 1 SUSY without a central extension will be identi- 
fied with the general eight vertex models with an external field. If the Boltzmann 
weights of the general eight vertex model satisfy a 'free fermion' condition, the 

model is exactly solvable and the free energy was derived first from dimer method 
[29] and later diagonalizing the transfer matrix [30]. Also, this model has been 
identified with the general XY-spin chain model with an magnetic field [31]. This 
relation with the lattice model will be very useful in our derivation of TBA equa- 
tions for the SShG model. It turns out that the SShG model is at the critical point 
of the XY-spin chain model. 

We organize this paper in the following order. In next section, we write down 
the lagrangian of the SShG and SSG models and derive the energy-momentum 
tensor supermultiplet and their relations under the SUSY transformation. Also 
we present the ^'-matrices of the models. In sect. 3, we review the basic formulae 
of TBA analysis which will be used in the next section where explicit derivations 
are explained. In sect. 5, we construct the FFs of the SShG model using the Wat- 
son equations and SUSY relations of the energy-momentum tensor. With exact 
two-point FFs, we derive the UV central charge of the model using the spectral 
representation of the C-theorem. 
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2. iV = 1 SUSY Integrable Model and Factorizable ^-Matrix 



We present the energy-momentum tensor supermultiplet of the N — 1 SSG 

and SShG model and the ^-matrix of the theories. 

2.1. Lagrangian and Energy-Momentum Tensor 
We start with a langrangian of a general N — 1 SUSY 



£($) = ^D<^>D<^ + iW{<^>) 



(2.1) 
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with a scalar superfield $ 



9) ^ (f) + iOTp + i\eeF, (2.2) 
and D and -D, the covariant derivatives 

Da = -=« + ^(7'^^)a5/.• (2.3) 
oQ 

The Grassman variable ^ is a Majorana spinor* In terms of the component fields, 
one gets 

m = l{d,<pf + \^[^ + W"{ct>)\ ^ + \ [W'{<t>)f. (2.4) 
The SSG model is a particular case of Eq.(2.4) with the superpotential 



777 

^($) = _cos(/3$). (2.5) 



The SShG model is the same superpotential with the purely imaginary coupling 
constant (5 = i(5. The A'" = 1 SUSY algebra is generated by the conserved charges 



-k Dirac matrices are 7 = 7 = 

' 1 oy Vi 
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Qi and Q2 

Ql^p+, Ql^P-, and {gi,g2} = o, 



(2.6) 



with the hght-cone momenta defined as P± — E ± P. These charges act on the 
component fields by 



Qi(f)^iijji, QiV'i = d+(f>, = F, 

(2.7) 

Q24> = #2, Q2^2 = -d-4>, Q2i'i = -F, 



with F = -W'{(f)). 

Integrabihty of the SSG and SShG models is estabihshed because they are 
equivatent to Toda theory based on the twisted super affine Lie algebra C(^)(2) 
[32-34]. The equations of motion of the SSG theory can be rewritten as super 
zero-curvature conditions. An infinite number of conserved charges at the classical 
level were derived [35] and checked to be preserved at the lowest order quantum 
corrections [36]. 

The energy-momentum tensor supermultiplet can be expressed by [37], 

Ja^= [(^$-W'(*))7/.^*]«, (2.8) 
or in light-cone coordinates 



with x± = h{xi ± xq) and d± = di ± Oq. In terms of the component currents, 
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one gets the energy-momentum tensor of the SSG model 



T+- = T_+ = sin2 /30 - !^ cos/30, 



(2.11) 



and its superpartner 

TTl 

^^1+ = ii>id+(j) ^'2+ = sin/30 

m (2-12) 

Including an appropriate normalization factor of 47r, we define the following nota- 
tion for the SUSY energy-momentum tensor: 



T = AttT. 



r = 47rr__, e = 47rr+_, 



(2.13) 



and their SUSY partners, 



Tf = 47r*i+, rF = 47r*2-, ©f = 47r*i_, ©f = 47r*2-h- 



(2.14) 



They are related to each other by the SUSY transformation 




(2.15) 



Q2TF^2iT Q2T^-\d-TF g2©F = 2i© g2© = - ^^-©f. 
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2.2. On-Shell Particle States and ^-Matrix 



If the coupling constant of the SSG model in Eq.(2.5) becomes pure imaginary, 
we have a simplest N = 1 SUSY field theory, namely the SShG model. Since the 
potential is not periodic, the soliton spectrum does not exist any more and the 
spectrum consists of only the fundamental particles appearing in the lagrangian, 
one scalar and fermion supermultiplet. We will denote on-shell states of these 
particles by \b{9)) and |/(^)) with a rapidity 9 which is related to the momentum 
by = m cosh 9 and P — m sinh 9. 

The SUSY charges defined in Eq. (2.6) can act on on-shell states as (See sect. 5. 3) 



me 



'/'\b{9)), Qr\b{9)) 



me 



e/2 



|/W>, 



Qi\fiO))- 

Q2\f{9)) = -zV^e-'/^\b{9)), Q2\b{9)) = zV^e-'/'\f{9)). 



(2.16) 



It is easy to see that this satisfies N — 1 SUSY algebra, Eq.(2.6). Action of SUSY 
charges on multiparticle on-shell states can be easily worked out using this and the 
anticommutivity of Qa and the fermion. 

Exact ^'-matrix of the SShG model was derived using the Yang-Baxter equa- 
tion, unitarity and crossing symmetry along with the commutativity of the SUSY 
charges and the -S-matrix [28]. In the basis of two-particle on-shell states in 
the order of \b1b2) , \f1f2) , \b1f2) , \f1b2)* the ^'-matrix has been obtained to be 
{9 = 91-92): 



S{9) = Y{9) 



V 



2i sin QTT 
sinh 6* 

i sin an 
cosh I 






cosh ■ 



2i sin an 
sinh 9 









1 








sinhf 



(2.17) 



sinh ; 



1 / 



with an arbitrary constant a which will be related to the coupling constant /3 of 
the SSG model in a moment. The prefactor Y{9) is needed to make the S matrix 



★ We use a short notation \b1b2) = \b{6i)b{62)) etc. 
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unitary and crossing symmetric. The following integral form will be useful later: 



Y{e) 



sinh I 



sinh I + i sin(|Q;|7r) 



exp — 



dtsmh{\a\t) sinh((l - \a\)t) ^.^^ 9t 
t cosh^ I cosh t TTi 



(2.18) 



With Y{9) = Yii-K — 9) and a factor of i arising in the crossing relation for hh — > // 
channel, the 5*- matrix of Eq.(2.17) is crossing symmetric. 

To determine the constant a we should refer to another derivation of the SSG 
breather jS-matrix. Using the SSG soliton jS-matrix, one can compute four soliton 
scattering amplitudes. By taking bound state poles of the incoming and outgoing 
soliton-antisoliton pairs one can derive the ^'-matrices of the SSG breathers [27]. In 
particular the ^'-matrix of the lightest breathers is the -S-matrix of the fundamental 
particles of the SSG and SShG models: 



where 



S{9) = Y{9) ■ n{9) ■ Sq{9) 
sinh 9 + i sin(2Q;7r) 



So{0) 



n{9) 



sinh 9 — i sin(2a7r) ' 

/ 2 _j_ 2i sin OTT 



sinh^ 

sin an 
cosh I 






cosh I 



-1 + 



2i sin Q7r 
sinh 6* 











1 








sinh I 



(2.19) 



sinh : 



1 / 



The factor 5*0 is the lightest breather ^'-matrix of the SG model. The constant a 
in Eq.(2.19) is given by the coupling constant of the SSG model [27] 



7 

167r 



-/32/47r' 



(2.20) 



For the SShG model with (3 — i(3 {(3 real), this constant reduces to 



a — 



l + ^/47r' 



(2.21) 



and — 5 < a < 0. 



11 



Two S'-matrices, Eqs.(2.17) and (2.19) arc equivalent. The sign difference in 
tfie // ff channel is explained l:)ecause all particles are considered as bosons in 
Eq.(2.19) by including the exchange factor —1 arising in // ff in the S*- matrix 
element. In this convention, the crossing relation is satisfied without any extra 
factor because all particles are bosonic. Besides, for the SShG model with a < 0, 
the 5*0 has no pole in the physical strip. Therefore, -S'o is nothing but a CDD 
factor and can be removed under the minimality assumption. For the SSG model, 
however, with a coupling in < o; < ^ < 47r/3) the Sq does have a bound 
state pole corresponding to the second breather. 

For a complete description of the SSG model, one should include all the S'- 
matrices of the solitons and breathers as was done in [27]. Depending on the 
values of the couphng constant of the SSG model, the spectrum of the bound states 
changes. In particular, if the coupling constant is in the range of ^ < ^ < 1, only 
the lightest bound states can exist along with the soliton and antisoliton in the 
spectrum. If the solitons are truncated from the theory keeping only the lightest 
bound states, the scattering theory becomes perturbed CFT by the least relevent 
operator. The UV CFT is the SUSY extension of the Yang-Lee model [26,27]. The 
^'-matrix is given by Eq.(2.19). 
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3. Thermodynamic Bethe Ansatz 



3.1. Diagonal TBA 

The TBA computes the Casimir energy for a theory on a circle of length R 
with 5"- matrices and particle spectrum as input data [38-42]. With a temperature 
T = 1/R the configuration of minimizing free energy gives the ground state energy 
of the system, which is again related to the central charge of the underlying UV 
CFT by 

E{R) ~ (C - 12Amin - 12A„,in) , (3.1) 

as R ^ (or T oo). Amin (and Amin) stands for the lowest conformal dimension 
allowed by the conformal theories. For unitary theories Amin is zero for the identity 
operator while it is in general negative for nonunitary theories. 

Consider 1 particles in a box of length L with periodic boundary condition 
(PEG). If we move the (A^ + l)-st particle of mass rria and rapidity 6;^ exchanging 
with all the other particles and come back to the original configuration, we get the 
following PEG equation: 

N 

^im.Lsu.l.e,Y[SaaM-Oi)^l, (3.2) 
1=1 

where the index specifies species of the i-th particle. For the diagonal scattering 
theories, the product of ^'-matrices are just C-number and one can take logarithms 
on both sides to get 



rUaL sinh + ^ - In Saa^ {Ok - Oi) = 27mk, (3.3) 



i=l 



with an arbitrary integer rik- These transcedental equations give solutions for 
{6i} for a given set of arbitrary integers {n^}. Therefore, considering all possible 
integers, and ets N, L — > oo the solutions form a band structure and one can 
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introduce a density of the rapidity states, Paid), defined by the number of allowed 
rapidity states between 6 and 9 + dO divided by dO. Any N rapidities of these 
states can be a solution of the PBC. 

Therefore, Eq.(3.3) can be expressed by 

27rpa(^) = maLcosh^ + I dd' pl{e')<Pab{e - e'), (3.4) 

b 

where the Pa{0) is the density of rapidity states which are actually occupied by the 
on-shell particles and 

<Pab{0) = -^g^lnSab{e). 

Introducing 'psuedo-energy' Ca defined by 



Pa{9) l + e-^'-W 
one can express the ground state energy by 



(3.5) 



E{R)^-Y,maJ ^cosh^ln(^l + e-^''(^)) , (3.6) 

-oo 

where is determined by the minimizing condition of the free energy: 

00 

eaie)^maR cosh e-J2 I ^<PabiO-0')^n(l + e-''^^'^). (3.7) 

b -oo 

This TBA equation can be solved easily for the UV {R 0) and IR {R oo) cases 
because the equations can be effectively described by simple algebraic equations. 
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3.2. NONDIAGONAL TBA 



For the nondiagonal theories, the product of ^'-matrices in Eq.(3.2) is the 

monodromy matrix and the PBC equation can be expressed as 

e'"^"^^^'^^^r„«(^|^i,...,^A,) = l, where 

Tabm. • • • , ^iv){::} = E sx'^{e - e,)szi{e - ^2) • • • - ^iv)- ^^-^^ 

If we add these equations for the index a, we can express it in terms of the transfer 
matrix 

^imL^mhej,^Q^Q^^ . . . , ^A^) = Nc, (3.9) 

where the integer Nc is the number of colors and the transfer matrix T = ^ Taa 
acts on V®^ . Precisely speaking, this is 'inhomogeneous' transfer matrix because 
it depends on each rapidity of in-coming particle states. 

To derive the TBA equations, one must diagonalize the transfer matrix which 
is quite a difficult task due to the size of the matrix and the inhomogeneity. In 
the lattice models, many pretty ideas have been invented for the purpose [14]. 
Although these methods are, in principle, applicable only to homogeneous cases, 
some minor corrections make it possible to apply it to the inhomogeneous ones. 
Among these, two methods have been successfully applied to derive TBA equations. 
The first one is using the inversion relation of the transfer matrix. TBA equation 
for the RSOS model corresponding to N — 1 SUSY sohton scattering theory has 
been obtained in this way [8] . Due to the periodic property of the transfer matrix 
and its inverted matrix, the eigenvalues are completely fixed by the location of 
'zeroes'. These zeroes satisfy constraint equations in terms of the rapidities ^^'s. 

Another useful approach is based on the algebraic Bethe ansatz method. In 
this method, one construct the eigenstates in terms of the monodromy matrix 
element acting on the vacuum. The eigenvalues can be directly obtained with an 
additional contraint which makes the eigenstate ansatz to be actual eigenstate. For 
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the example, this method can be used for the SG model and N = 2 SUSY theories 
[21-23]. 

In general, the eigenvalues of the transfer matrix have the form like 

N 

Aie\ei,...,eN)^l[fie-xk), (3.10) 

k=l 

with the zeroes of a function f{x) (/(O) = 0) which satisfy 

N 

l[9{xk-ei)= Const. (3.11) 
1=1 

By taking an imaginary part of the logarithms of both sides of Eq.(3.8) and intro- 
ducing pscudo-cncrgics both for the real particle states and for the zeroes, one can 
find the TBA equations which are very similar with those of the diagonal TEA. 
Only difference is there is no mass term ('driving term') in the TBA equation for 
the density of the zeroes. This absence of the driving term makes a big difference 
in the analysis of the TBA equations. 

3.3. Casimir Energy and UV Central Charge 

The TBA equations in the UV limit {R 0) can be easily solved because the 
pscudo-cncrgics become independent of the rapidity around = 0. This plateaux 
extends upto 6 ~ — ln(mi?). Therefore, the pscudo-cncrgics become practically 
constant in the limit and the TBA equations can be reduced to the mere algebraic 
equations like 

Xa = ll{l+Xbf^\ (3.12) 
b 

with Xa — exp(— ea(0)) and 

Nab = ^ [0a6(oo) - <f>ab{-Oo)] . (3.13) 

If all mass terms are non-zero, the €^(00) diverges like ea{9 00) ~ uiaR cosh 9 
and Ua — exp(— 60(00)) vanishes. If some of the driving terms are zero as is the 
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case for nondiagonal TBA, ya' satisfies 



b' 



(3.14) 



where a' denotes species with vanishing driving term. The ground state energy can 
be expressed compactly with these variables Xa and Ua by 



Xr, 



-£ 



Va 



1 + ya 



(3.15) 



where J0,{x) is the Rogers dilogarithmic function 







ln(l-i)_^ \nt 



(1-t) 



(3.16) 



From Eq.(3.1) the central charge of the UV CFT is given by 



C - 12(A„iin + Amin) = ^ 



Xn 



(3.17) 
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4. TBA for the iV = 



1 SUSY models 



In this section we diagonalize the inhomogeous transfer matrix associated with 
the SShG ^'-matrix. The essential observation is that Eq.(2.19) satisfies so-called 
'free fermion' condition of the general eight vertex model with an external field. 
We will derive the TBA equation based on the inversion relation for the transfer 
matrix. We apply this TBA equation to both the SShG model and a SYL model 
perturbed by the least relevent operator and derive correct UV central charges. 

4.1. Free Fermion Models 

After the celebrating solution of the symmetric eight vertex model by Baxter, 
Fan and Wu obtained an exact expression of the free energy for the general eight 
vertex model with an external field if the Boltzmann weights satisfy some additional 
contraint, named the free fermion condition [29]. They called this model 'Free 
Fermion' model (FFM) although the name is slightly misleading. The model turned 
out to be highly non-trivial and interacting. 

We start with the Boltzman weights of the general eight vertex model: 

/a+ d \ 

6+ c 

c 6_ 

\ d a-/ 



(4.1) 



for the following vertex configurations: 



T 

a+ 



T 
d 



_T_ 
i 

d 



(4.2) 



If R{9) satisfies the Yang-Baxter equation and the free fermion condition 



a-i-O- -|- b+b- — c + d , 



(4.3) 



and if the following combination of the Boltzman weights are independent of the 
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rapidity 



^ - o.. . , . \ A ^ (4-4) 



a+b- + a-b+^ 2{a+b- + a-b+) ' 

the transfer matrix T commutues; [T{u),T{v)] — 0. Due to this commutativ- 
ity, there exist infinite number of conserved charges including a Hamiltonian of 
the corresponding one-dimensional spin-chain model. This Hamiltonian has been 
identified with that of the XF-model with a magnetic field, 

N 

HxY = + + ri^'^Z+i + ^/'^Z+i) - K]' (4-5) 

where cr^ = ^(cr^ ± ia^) with a conventional Pauli matrices. 

To identify the FFM with = 1 SShG model, we rewrite the jS-matrix of 
the SShG model, Eq.(2.19), by rearranging the two-particle basis. In the order of 
|//), the i?-matrix of Eq.(2.19) becomes the general form of the 

FFM with 

2i sin ttTT i sin an , sin ott , ^ , 

a± = ±H r-r-7-> &± = 1, c=— -5-, d= —q, 4.6 

smh 6' smh | cosh | 

if we identify t and ^ with |6) and J, and with |/).* 

It is an easy exercise to check these weights satisfy the free fermion condition 
Eq.(4.3). Also, the constants F and h becomes 

r = sinQ;7r, h = —1. (4.7) 

Since /?. = — 1 is a critical point of the XF-model, the SShG model corresponds to 
the critical point of the general eight vertex model with free fermion condition and 
with vanishing modulus. 



★ The Boltzman weights in Eq.(4.2) become the S'- matrix element if we adopt a convention 
that time flows from bottom-left to top-right (/'). 
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4.2. Inversion Relation 



Critical step to derive TBA equations for the SShG model is to find an inversion 
relation for the transfer matrix. We derive the following inversion relation in the 
Appendix A: 

T{u\9i, ...,dN)T{u + ni\di,..., 9n) = (-l)^x 

■ N N / N N 

i[M+iu - +\[m_{u -ei)+F ( n^+(^ - + n^-(^ - 

.1=1 i=l \i=l i=l 

(4.8) 

where the fermion index operator F is either +1 for the bosonic state or —1 for 
the fermionic one. 

The functions appearing in Eq.(4.8) are expressed in terms of the Boltzmann 
weights as follows: 

M-i- = a+a_ — d , M_ = a+a_ — c , 

= sinh^ (j) a+b+ + cosh^ a_6_ — 2 sinh (p cosh (f) cd, (4.9) 
F_ = — cosh (j) a_|_6_|_ — sinh a_6_ + 2 sinh cosh cd, 



and 



2cd 

tanh(2(/)) = — — = sin an. (4.10) 

^ ^ a+b+ + a-b- ^ ' 



Using Eq.(4.6) one can find 



sinh (I + ian) sinh (| — ian) cosh (| + ian) cosh (| — ictTr) 

= —a —a , — —a — 3 , 

smh I smh | cosh 2 cosh 2 

^ cosh (5 + ia'n) sinh (| — ia'n) ^ sinh (| + ictTr) cosh (| — iojTr) 
cosh I sinh | ' sinh | cosh | 

(4.11) 

From these expressions one can notice that under the change u ^ u + in 



M± and F± F^, 
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therefore, T{u)T{u + Tii) = T{u + Tii)T{u + 27ii). This means 



(4.12) 



These matrix relations can be easily transformed to equations of the eigenvalues 
of the transfer matrices; A{u\9i, . . . , 9n) is 27ri symmetric function, 



A{u\9i, ...,9n) = A{u + 27ri\9i, ...,9n), 



(4.13) 



and the inversion relation is nicely factorized. 



A{u\9i, 9n)A{u + 7ri\9i, ...,9n) 

-itcosh(^ + z|a|7r) -iX- sinh + i|a|7r) 
11 ^^oV. fu-9i\ +^11 



=1 cosh (^) 



l\ sinh (^) 



X 



4= 

iV 



^ cosh (^^ - i\a\n) -A- sinh (^^ - i\a\TT) 



n 



J cosh (^) 



J sinh (^) 



(4.14) 



4.3. Eigenvalues 

Since A(m) is a 27ri-pcriodic function with poles at m = 9k and u = 9k + ill, it 
can be completely fixed by the location of zeroes on the strip in the complex plain 
of —in < lm[u\ < in and — oo < Re[u\ < oo. Also from the fact that A{u) becomes 
a constant as — > oo, we can find that 



N 



A{u\9i,...,9n) = Const. JJ 



sinh 



sinh 



sinh (^^) cosh (^^) 



(4.15) 



where the 2A'" zeroes {z'^} and {zj^ } located on the strip will be determined as 
functions of ^^'s. 
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We defined the zeroes in the way that z'^ and z^ come from the first and 
second factors of the RHS of Eq.(4.14), respectively. Therefore, they satisfy 

N tanh f ^tfi + N tanh f - iblTr) 

n — ^ ,,J. n — ^ 

The solutions of these equations can be written in terms of real variables in the 
following way: 



z'^ — Xk — i\oi\'K^ — i\a\'K -\- 

2;^ = Xjfc + i|Q;|7r, xif + i\a\T{ — 

where a real number xj^ satisfies 



(4.17) 



N tanh 



Out of four possible choices of z^j: for k = 1, . . . ,N from Eq.(4.17), only two 
choices arc allowed. This can be understood easily if one considers the limit of 
\a\ 0. The Boltzmann weights arc cither +1 or —1 from Eq.(4.6). This means 
the transfer matrix is just a constant matrix without any dependence on the ra- 
pidities. Now from Eqs.(4.15) and (4.18), the only possibility for the eigenvalues 
of the transfer matrix to be independent of 9i is when {xj^} — {9i} and when 
^k~^k ~ -'-^^ ^- example, if one choose {z^, z^) — {9]f—i\a\TT, 9k+i\a\TT) 

pair for some k as \a\ — > 0, the eigenvalue will get term like tanh (^^^-^^- Obvi- 
ously, this eigenvalue should be excluded for the constant transfer matrix. This 
leaves only two choices for the zeroes: 

i^k ^ ^k) ~ (•'^k ~ 'i\(^\''^j^k + ~ i''^) {xjf — i\a\7T + in,xi( + i\a\7i). (4.19) 



From the product form of Eq.(4.14), one notices that if we choose one pair of 
zeroes in Eq.(4.19) the other pair become zeroes of A{u + i7r). Since one can choose 
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in - i\a\Tr 

o«o«on«o«o« 



> 

• 0«0«(>0«0«0 

— i|Q;|7r 

Figure 1. The zeroes on the complex 9 plane. 

the zeroes between the two possibihties for each k {k = 1, . . . , N), we can construct 
2^ different eigenvalues in this way. Also, one can convince that Eq.(4.15) satisfies 
Eq.(4.14) because if we dividing the RHS of Eq.(4.14) with A{u)A{u + iir) using 
Eq.(4.15) the final expression has no poles and zeroes while it is bounded. This 
means the ratio should be a constant. 

Using all these results, the eigenvalues are compactly expressed by 
A(^i)ei,...,€Ar = Const. 

with 

X,{9) = smh + 6^ j cosh i^- - e^j . (4.21) 

€ — +1 corresponds to the first choice in Eq.(4.19) and e — —1 the second. (See 
Fig.l) The real zeroes Xk are determined by Eq.(4.18). 



Hill sinh(M - 6i 



= ±1, 



(4.20) 
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4.4. NONDIAGONAL TBA 



From Eq.(3.8) and Eq.(4.20), the PBC equation becomes 



im sinh 9 



N 

n 

i=l 



Y{9 - 9i) 
smh{e - Oi) 



N 



\{K,{e\^U...,XN) = l, (4.22) 
k=l 



and the constraint equation (4.18) in the hmit N,L ^ oo are expressed by intro- 
ducing the densities p for the allowed states, for the occupied states, P_|_ for 
e — +l zero-state, and P_ for e = — 1. In terms of these densities, one gets 

27rp{9) — m cosh 9 

+ J de'[p\e')^Y(e - e') + p+{e')^+{9 - e') + p_(^')^-(^ - o% ^^^s) 

2nP{9)^ J d9'p^{9')^T{0-0'), 



where 



d_ 



Im In 



1 d 



y{e) 

sinh 9 



tanh ( I 



^±(^) = ^InilnA±(^), 



tanh (I + ^) 



(4.24) 



Using Eq.(4.21) and A_ = (A-|-)*, one can easily show that the kernels are 
related by 



$t(^) = 2$+(^) = -2$_(^) = -4 In 

I uu 



sinh 9 — i sin 


a 


71 


sinh 9 + i sin 


a 


7r_ 



(4.25) 



which is nothing but the kernel of the sinh-Gordon model. We will denote this 
kernel by Also, we can eliminate P- from the first equation of (4.23) using the 
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second equation and P = P-f + P_ to rewrite it as 



2Tip{9) =mcosh9 + J dO' 



p^{e')[^Y - * $](^ - e') + p+{e')^{e - e') 



.26) 



where we introduce a convolution defined by 



dd' 

2^ 



f{9-e')g{9'). 



Eqs.(4.23) and (4.26) have the same form as those met in diagonal TBA hke 
Eq.(3.7), except that the second equation in Eq.(4.23) has no driving term. The 
TBA equations, therefore, can be written down as before, 

mR cosh = e(^) + - 1$ * $] * ln[l + e"']^ (0) + * ln[l + e"'^]) (0) 

= £(^) + ($*ln[l + e-^]) (0), 

(4.27) 

in terms of the pseudoenergies e and £ defined by 

p\0) e-<^) P+{0) e-^W 



p(0) i + e-^W P{0) l + e-^W 



(4.28) 



4.5. Central Charges of the SShG model 

In the UV limit, Eq.(4.27) reduces to simple algebraic equations of the variables 
X = exp[— e(0)], X = exp[— £(0)] as argued before. For the SShG model, the 
algebraic equations become 

x^{l + xf{l + X)\ X^{l + x)\ (4.29) 

with 

oo oo 
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It is not difficult to compute these exponents using Eqs.(4.21) and (4.24), 



oo oo 

f dd 1 f dd 



and 



oo 

/ 



2^ 



-1 2 



d9_ 
2^ 



m 



a = and 6 = 1. 



(4.31) 



Using these values, the solution of Eq.(4.29) can be found easily as 



X — oo and X — oo. 



(4.32) 



One also needs the psuedo-energies as ^ — > oo. Since the mass term for e(^) is 
non-zero, e diverges as ^ — > oo, and y — exp[— e(oo)] = 0. Then, from the second 
equation of (4.27), Y = exp[— £^(oo)] = 1. 

Using all these values in Eq.(3.17) the UV central charge becomes 



6 



£(l)+£(l)-£ 



3 
2' 



with Ajnin = Amin = for the SShG model. This is correct UV central charge of 
the SShG model with a boson and a fermion. 
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4.6. Central Charges of the SUSY Yang-Lee Model 



As explained in the previous section, one can truncate all solitons from the 
SSG multi-sohton Hilbert space to have only breathers. In particular, for the 
coupling constant a = ^, only the lightest breather and its supcrpartner can exist 
in the spectrum with the .S-matrix given in Eq.(2.19) [27]. This is the SYL model 
perturbed by the least relevent operator. 

The fundamental difference from the SShG ^'-matrix is that because of a > 
So{9) is no more CDD factor. It has a pole which is identified with the particle 
itself. If we denote the particles as B and F, the bootstrap relations are 

BB (FF) ^ B^BB {FF), BF (FB) ^ F ^ BF (FB). 

Except this difference, all the TBA analysis of the SShG model is equally applicable 
to the SYL model. 

The SYL conformal theory can be constructed as a coset CFT given by 



The central charge of the model is C = Due to the nonunitarity of the model, 

the lowest conformal dimension is not zero. To determine the lowest conformal 
dimension of the model, we refer to the general formula of the general coset theories. 
The primary fields of the coset theory have the conformal dimensions given by the 
following formula [43] : 



with the restrictions, 

0<1<K, l<r<p-l, l<s<p'-l, I ^ \r - s mod 2K\. 

For the SYL model with the values oi K — 2, p — 2, p' — S, one finds that 
only r = 1, Z = 0, 1, 2, and s — 1, 2, 3, 5, 6, 7 are allowed. The minimal conformal 



SU{2)k^SU{2)l 
SU{2)k+l 



K^2 and L + 2 = -. 



_ l{l + 2) f {^rp'-spf-{p'-pf 
A{K + 2) 4K AKpp' 




p' ^p + Kq, 
(4.33) 
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1 

'4- 



dimension arises when (8 — 2s)^ is minized, i.e. with s = 3. Therefore, the minimal 
conformal dimension of the SYL model is Amin = 3 

Now we compute the Casimir energy of the model using TBA. Notice that the 
only change from the SShG TBA is that the kernel $y in Eq.(4.27) gets extra 
factor — ^ln5'o(^) due to 5*0 in Eq.(2.19). With a > this introduces extra — 1 
in the exponent a in Eq.(4.29) to make a — —1. With this change the algebraic 

equations now become 

l+X 



X 



X ^1 + x, 



(4.34) 



and the solutions are x — -\/2 and X — 1 + \f2. Inserting this into Eq.(3.17), one 
can get 



6 



7r^ 



^/2 



1 + ^2 



'1W2' 
2 + ^2 



-L 



3 

4' 



which is exactly C — 24Aniin as required. 
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5. Form Factors of the Supersymmetric Theories 



In this section, we derive FFs of the SShG model from the 5'-matrix of the 
modeL Starting with the axioms for the FFs to satisfy, we write two-point FFs in 
terms of two factors. The first factor is determined by the properties of operators 
and contents of poles while the second factor is completely determined by the 
eigenvalues of the ^'-matrix. For the SShG model, we compute one-point FFs first 
to fix overall normahzation. Then, two-point FFs will be derived using the Watson 
equations and SUSY relations of the FFs. 

5.1. Form Factor Axioms 

We start with axioms of the FFs [6]. Denoting \a{e)) as an on-shell particle 
state of type a with a rapidity 9, a matrix element of an Hermitian operator O 
between vacuum and in-coming states can be expressed by 

(^1, ■■■,0n) = (0|O(0)|ai(^i), . . . , an(^n)>i„ (5.1) 

with a normal ordering of rapidities, 9i > 02 > ■ ■ ■ > On- This FF should satisfy 
the following axioms: 

1. F is analytic in each variable 9ij = 6i — 6j inside the strip of < Im[^y] < 27r 
except for simple poles. All other type of FFs can be reduced to the standard form 
hke 

■ ■ ■ , bm{e'J\O{O)\ai{0i), . . . , an{0n)\^ 

with the charge conjugation matrix C arising in the crossing process. 

2. Due to the normal ordering, any exchange of two rapidities should involve 
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a scattering process; 



-^ai,...,ai,ai+i,...,a,i(^l) • • • 5 ^15 ^i+l) • • • ) ^n) 

. . ' ' 5.3 

with the phase factor arising from the exchange of particles {fa = for a boson or 
1 for a fermion)* 

3. Relativistic invariance dictates 



Fa,,...,ar.{01 + A, . . . , + A) = e'^Fa,,...,aM, • • • , ^n), (5.4) 

where s is the spin of the operator O. 

4. A 27ri translation of one of rapidity is equivalent to 

Fa,,...,ar.{ei,...,en + 2m) = (-1)^- ^ Fa„,ai,...,a„_i (^n, ^1, • • • , ^n-l) • (5.5) 

There are two origins of the poles. One is the annihilation pole and the other 
corresponds to the bound states. Existence of the poles give extra constraints on 
the FFs. 

5. The annihilation pole arises in the channel of two incoming particles which 
are related to each other by the charge conjugation at ^ = ztt. This gives 

Res 

^0' ^ QFa^,...,a„,a,a{^l, ...,0n,0' + in, 0) 

= [1 - Taam, On)t':::X ^«i,...k(^i, • • • , 

using the monodromy matrix defined in Eq.(3.8). This equation relates a n + 2- 
point to a n-point function. 



★ We do not need this phase factor if we use Eq.(2.19) as the S'-matrix instead of Eq.(2.17) 
since all particles are bosonic. 
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6. The bound state pole of a S'-matrix with a residue 



i Si^\e) = -Vl,vf\ (5.7) 



the FF should satisfy 



Res 

a n ^ai,...,an,aa,ab{^l, ■ ■ ■ ,0n,0a,db) — ^ab^ai,-,an,ac{^l, ■ ■ ■ ,0n,dc)- (5i 



5.2. Two-point Form Factors 

The general axioms can be much simplified for two-point FFs. It can be written 
in terms of two factors like, 



-P'oi02(^l) ^2) = -f^ai02(^l) ^2)-^^a'^(^l ~ ^2), 



(5.9) 



where F™^" satisfies the Watson equation without any pole, 



(5.10) 



and the prefactor i^^ai02(^i) ^2) has all the required poles and operator dependence. 
Note we omitted a phase factor from the fermion exchange operator treating all 
particles as bosonic. 

can be determined from the following steps: In the basis which diagonal- 
izes the S'-matrix, Eq.(5.10) becomes a simple functional relation. Then, using an 
integral representation for the i-th eigenvalue of the ^'-matrix, one finds 



Si{9) = exp 



dt ^ , ^ . Bt 
— /^(^)smh — 
^ m 



exp 



dt h{t) . 2 9t 

r-TTT sm — 

t smh.9 ZTT 



LO 



, (5.11) 



where 9 — in — 9. Rotating back to the original on-shell two-particle states, one 
finds the F™^. 
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The function K{6i,62) should be determined by the other axioms. Eq.(5.5) 
requires it be a symmetric function of the rapidities if ai = 02 and has a ivr-pole 
if ai = 02. The asymptotic behaviour under the rapidity translation is related to 
the spin of the operator by Eq.(5.4). The overall normalization of the FFs is fixed 
by the one-point function. 

5.3. One-Point Form Factor of the SShG model 

We work out one-point function of the SShG model to fix normalization of the 
general FFs. From the Fourier transformation of the elementary fields 

dk 1 



0(x) 



2vr y^ko 
dk 1 



(5.12) 



with u* — V and with the commuation relations 



{fkjl]=27rko6{k-k'), 



27rko6{k-k'). 



(5.13) 



On-shell SUSY is determined from the SUSY transformation of the elementary 
fields, Eq.(2.7) and \b{e)) = b^e)\0) and |/(^)> = fH0)\O): 

■ P+ 



Qi\m) 



■\b{9)) Qi\bi9)) = tV2vi{9)\f{9)) 



V2vi{9) 

Q2\f{e)) = ^^j^^\m) Q2\h{9)) = iV2v2{9)\f{9)). 

From = (0|ga [0(0) |/(^)>] , one gets 

Qa\f{e)) = -i^ua{9)\b{9)), 
and comparing this with Eq.(5.14), one can find the spinors 



(5.14) 



(5.15) 



This gives the SUSY transformation of on-shell states given above in Eq.(2.16). 
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Combining this and Eq.(5.12), we fix the one-point function as follows: 



2 \ 



(5.16) 



5.4. Two-point form factors of the SShG model 

We compute two-point FFs of the trace of energy-momentum tensors and their 
SUSY counterparts, and ©p (©f)) given in Eqs.(2.13),(2.14). These operators 
are of particular interest for their role in the C-theorem. First, we derive SUSY 
relations between the FFs using Eq.(2.15). 
From = (o|(5a[(^|ai(^i)a2(^2))], one finds a relation 

<0|ga[O]|ai(^i)a2(^2)> - -(-l)^(^\0|C»|ga[ai(^i)a2(^2)]>, (5.17) 



with F is 1 for fermionic and for bosonic. This gives the following relations 
between the FFs: 



hb-^ 2 










1 


Fbb- 2 







hi 



Of 



bf 



2 

m 



~iF, 



Of 
bf 



X^F 



fb 



-1 



pOf I 
^bf + 



-F 



Ob 



fh 



(5.18) 



where each FF is a function of Oi and 6*2 or of Xi = e \ 

A special case of a — > 

It is useful to consider a case of a ^ where the S'-matrix is of diagonalized 
form of (1, —1, 1, 1) from Eq.(2.19). Let's compute F^ and In terms of the 
solution of the Watson equation* 



771 min 
^bb 



9 



Ff/- = sinh-, 



(5.19) 



7k- This solution is not unique in the sense that one can multiply any even function of 6 
satisfying /{in + 6) = /{in — 9). If we include these functions in the prefactor K, one can 
define F™'" uniquely. 
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the FFs can be written as 



Q 

F^{xi, X2) = Kj,i,{xi, X2), and Fff{xi,X2) = Kff{xi,X2) sinh -. (5.20) 



Since K should have the in pole (or at xi — —X2) and K{9i+A, 62+ A.) = K{9i, 62) 
because the spin of © is zero, we can find 



F^{xi, X2) = 27rm^, and F^(a;i, X2) = 27rm^ sinh (5.21) 

Here we fixed the normalization factor as nm? by comparing with the perturbative 
computation using Eqs. 

After finding these, one can derive the other FFs simply using Eq.(5.18) as 
follows: 



e 



F^/{xi,X2) = 27rm3/2y^, F^/{xi,X2) = 27rm3/2^. (5.22) 



One can check that the spins of Op and Op can be found correctly as under 
the rapidity translation. Also one can check that these FFs are consistent with 
pertrubative computation. 

For General a 

For general cases, we should diagonalize the S'-matrix first. Whether we use 
Eq.(2.17) with the phase factor in Eq.(5.3) or Eq.(2.19) without such factor, we 
find the eigenvalues of F = sector {bb and //) are complicated and hard to find 
the integral representations. Instead, we consider F — —1 {bf and fb) sector first. 
The S'-matrix is easily diagonahzed by the eigenvectors 

1+) = ^ (l^i/2> + |/i&2» , |-> = ^ (|6i/2> - |/l62» , (5.23) 
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with eigenvalues 



S+{9) = exp 



oo 

f dt ^ , , . Ot 
J t m 



S_{9) = - exp 



oo 

f dt ^ , , . Ot 
/ —f-(t)smh — 



LO 



(5.24) 



with 



^^^^^ ^ (l-c-oshf)(l + (-osh((l-2|n|)t)) ^ cosh((l - 2|a|)f) 



sinh^ t 



coshi 



(5.25) 



Prom these integrals, F^^^ in the basis of |+) and |— ) can be obtained as 



= exp 



dtf^it) . 2 Ot 
t smht ZTT 



LO 



, F^'^ = cosh - exp 



ot 

/ 



dtf-{t) . 2 ot 

sm — 

t smht 277 



LO 



(5.26) 



where we chose a normalization such that 



F± ^ 1 as a^O. (5.27) 

For numerical computations, we list expressions of F™™ which converge fast in the 
Appendix B. 

Now we consider FFs of Qp in the following form: 

F^^xi,X2) = K+{xi,X2)Ff'', F?^{xi,X2) = /^-(a;l,a;2)Fr^ (5.28) 

and similarly for Op in terms of K±. These K± and K± can be determined 
from the spins of the operators and symmetric properties of the states under the 
exchange, |+) — >• |+) and | — ) — >• — | — ) under xi X2, as follows: 

i^± = ^(V^±V^), K±^b(^±^], (5.29) 

where constants A, B can be determined by taking a ^ limit and comparing 
with Eq.(5.22). 
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Now rotating back to the on-shell states one can find 



F^/{x-L,X2)^27rm^/'' 



2^1^^ — ^ + V^— — ^ 



(5.30) 



Also from Eq.(5.18) one can obtain other FFs 



F§^{xi,X2) = 27rm^ 



^ + ^ cosh - 

2 2 2 



Ff f{xi,X2) = 27rm^ ^ ^ ^ ~ sinh -. 



(5.3i; 



Wc checked these FFs using the first order perturbative computations. FFs for 
other component of energy-momentum tensor can be written down by just multi- 
plying P+/P- to the above FFs. 
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6. Spectral Representation of C-Theorem 



We compute the UV central charge of the SShG model using the two-point FFs 
computed in the previous section and the spectral representation of the C-Theorem. 
This provides a consistency check for the FFs and shows the fast convergence of 
the FFs expansions of correlation functions. 

6.1. Spectral Sum Rule 

The C-theorem, first introduced by A.B. Zamolodchikov, plays an important 
role in the study of off-critical models [44] . The C-function, describing a degree of 
freedom of the 2D models, connects smoothly two renormalization group (RG) fixed 
points as the length scale of the theory increases from UV limit to IR. For some 
specific models like the perturbed minimal CFTs by the least relevent operator 
with positive coefficient, the renormalization group (RG) fiow connects two RG 
fixed points corresponding to two adjacent minimal CFTs [45]. This RG fiow will 
end up at the massive point with C — 0. 

This theorem can be neatly expressed in the following integral of the two-point 
correlation function of the trace of energy-momentum tensor following Cardy [46] : 

oo 

\x\>e (6.1) 

Ci(/x, A) = -^^Im 

Expanding the correlation function in terms of intermediate on-shell states, the 
spectral density function Ci can be expressed in terms of the FFs by 

Ci(/.,A) = ^^|(0|e(0)|a>|V(g-p,), (6.2) 

^ a 

where pa is the energy-momentum vector of the multi-particle state a and the 
vector q is defined as q — {fj,, 0). 
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For the massive theory, the sum rule of AC effectively gives the UV central 
charge since Cir vanishes. Although one needs the infinite number of the FFs to 
compute it rigorously, there are many evidences that the sum in Eq.(6.2) converges 
very fast for the massive theories [7-10]. With this observation, one can compute 
the UV central charge using the two-point FFs of © quite accurately. In next stage, 
we will compute this numerically using the FFs of the SShG model derived in the 
previous section. 

6.2. Sum Rule for the SShG model 

The two-point contribution to the sum rule becomes 



X S{m cosh 9i + m cosh $2 — ^)5{m sinh 9i + m sinh $2) j-g 3-^ 

00 

3 f de 



4 



with the FFs given in Eq.(5.31). 

For the special case of a = where the SShG model becomes free with a boson 
and fermion, one can insert Eq.(5.21) into Eq.(6.3) and using 



00 



de 2 r ,^sinh2^ 1 

de- 



r de _2 r 

J cosh^e ~ 3' y ''^cosh^^ ~ 3' 


one can easily find C — \. 

For the generic value of a we integrate numerically using the regularized ex- 
pressions for the F± in Appendix B. Using these we list AC^^^ as for several values 
of the coupling constant in Table 1. This shows a good agreement with the UV 
central charge C = |. The convergence of the SShG model seems slow compared 
with the sinh-Gordon result [9]. This suggests in the SShG model one arrives 
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1- 


n 


Ar(2) 


1 

999 


0.001 


0.9993 


1 
199 


0.005 


0.9953 


1 

99 


0.01 


0.9902 


1 

49 


0.02 


0.9800 


3 

97 


0.03 


0.9697 


1 

19 


0.05 


0.9495 


1 
9 


0.1 


0.9093 



Table 1. The first two-particle form factor in the Sum Rule of AC = |AC. 

strong coupling region earlier than the sinh-Gordon model as one can see from the 
fact that the limit of the SShG coupling constant is ^ = 5 while ^ = 1 in the 
sinh-Gordon model. 
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7. Conclusion 



In this paper we obtained two results on the = 1 SUSY integrable models. 

The first one is the computation of the UV central charges from TBA method. 
The nondiagonal TBA of the SShG and SYL models has been rigorously derived 
from the essential observation that the A^ = 1 SUSY models can be identified with 
the eight vertex free fermion models. These TBA equations produced correct UV 
central charges. 

The second result is two-point FFs of the SShG model using the FF axioms. 
Here the difficulty arising from the nondiagonal scattering theories has been avoided 
from the SUSY relations of the FFs. The spectral representation of the C-theorem 
showed that two-point FFs can give good approximations in the infinite sum of the 
intermediate states even in nondiagonal theories. 

Our results suggest some interesting directions to proceed further. Actually, we 
notice that wider class of A^ = 1 scattering theories are belonging to the eight vertex 
FFMs which will be reported in separate pubhcation [47] . The relationship between 
these SUSY models and the eight vertex FFM may have some deep structure 
because the FFMs seem to have interesting hidden symmetries [48]. In particular, 
it has been noticed recently that the FFMs have a hidden quantum group symmetry 
[49]. It would be interesting to see how this quantum group symmetry will be 
related to the N — 1 supersymmetry in the trigonometric limit. 

In this paper, we could not say much on the general FFs of the theories. The 
solution of the FF bootstrap equations are very difficult and are limited to only a 
few simplest diagonal theories. We can reduce, however, the nondiagonal bootstrap 
equations to the level of diagonal theories by diagonalizing the inhomogeneous 
transfer matrix. It will need some more work to solve these reduced bootstrap 
equations completely. 
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APPENDIX A 

Inversion relation for the Free Fermion Model 

Wc follow Fcldcrhof to diagonalize the transfer matrix of the FFM [30]. We 
want to point out, first, the difference of our derivation from the lattice model 
computation. The first difference is that we want to diagonalize the inhomogeous 
transfer matrix. This difference often introduces much difficulty for the computa- 
tion. However, this difficulty can be avoided by the second difference, which is that 
we are working with the FFM at the critical point. With this advantage, we can 
derive the inverse matrix of the FFM transfer matrix and, furthermore, express it 
using the original transfer matrix with slight change in the rapidity u. 

It is convenient to reexpress the Boltzman weights Eq.(4.1) in terms of the 
(7-matrices, 

fA{e) B{e)\ 

R(9) = \ , A.l 

^ ^ \c{e) D{e) I ^ ' 



A = a^a^a + b^a a^, B = da^ + ca , 
C — 07^ + cicr", D — b-a^a~ + a-a~a^ . 



(A.2) 



Then, the transfer matrix becomes 



T(«|^i,...,^Ar) =Tr2 



N 



X[R{u-ei) 



.1=1 



(A.3) 



Now we define new transfer matrix T\ corresponding to new Boltzman weghts 
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defined by 



a\. — —b±, b\. — a±, — c, and — —d. 



(A.4) 



In the same way as before, one can express Ti by 



Ti{u\9i,...,9N) = Tr2 



N 



l[Ri{u-9i 



.1=1 



Ri{9) 



'Ai{9) Bi{9) 
Ci{9) Di{9) 



, (A.5) 



where 



Ai — —hj^a'^a + a-f-U u''", Bi — —da'^ + ca , 
Ci = ca~^ — da~ , Di = a-a^a~ — h-a~a^ . 



(A.6) 



One can check that these new Boltzman weights again satisfy the free fermion 
condition Eq.(4.3). 

Next step is to show that TTi oc 1. For this purpose, we multiply two matrices 



T{u)Ti{u)=Ti2 



N 



.1=1 



Tr2 



TV 



i=l 



Tr 



2(g)2 



N 



.1=1 



(A.7) 



Defining the 4x4 matrix Ri Ri^i as Si, one can find a similarity transformation 
S', — XiSiX~^ where S[ is of triangular form. The X and S' are given by 



X 



cosh0 — sinh0 
— sinh0 cosh0 



V 73-71 ^ 



/ M+ * * * \ 

F_(7^ * 

F+a^ * 
\ 



M_ / 
(A.8) 

where M±, F±, and (j) are given in Eqs.(4.9) and (4.10). We did not specify the 
unnecessary non- vanishing components (*). 

Most important observation is that tanh^ becomes just a constant for the 
N — 1 supersymmetric theory. This means one can make all the SiS in the trace 
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of triangular form by the same similarity transformation X. Therefore, TT\ — 
Tr4 Jl^^ and from Eq.(A.8) one can derive 

- iV iV 

(A.9) 



T{u)Ti{u) 



.1=1 



N 



N 



\i=l 1=1 

with F — Ylcrf is either 1 (bosonic) or —1 (fermionic). 

Now, consider a translation u ^ u + iir. Under this the Boltzman weights of 
the SShG model change 



a± 



-a=p, b± — > 6=p, c ^ d, and d 



-c. 



(A.IO) 



Again this satisfies the free fermion condition. Now, the transfer matrix with 
translated rapidity can be expressed in terms of R2{u — 9) — R{u -\-m — 9) by, 

N - ^ 

^R2{u-9i) , i?2= I 7 1 , (A.ii) 



T{u + m) = Tr2 



.1=1 





B2\ 








D2) 



where 



A2 — —a+a'^a~ + b+a~a'^, B2 — —ca'^ + da~ —, 
C2 — da'^ — ca~ , D2 — h-u'^a' 
From Eq.(A.6), one can notice that 



a-G o""*". 



(A.12) 



(A.13) 



^2 = -^i, S2 = -Ci, ,C2 = -Si, and L>2 = -Ai. 
Considering the it!-matrices as 2 x 2 matrices, the Ri and R2 are related by 

i?2 = -(7*i?i(7^, 

where cr*' is the usual Pauli spin matrix. This gives 

T{u + ^ttI^i, ...,9n) = {-lfTi{u\9i, ...,9n), 
and from Eqs.(A.9) and (A. 14), the inversion relation Eq.(4.8). 



(A. 14) 
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APPENDIX B 



Regularized Expression for the Form Factors 

For the numerical computation we can rewrite F^^^ in Eq.(5.20) as follows: 





n 




oo 


F^^{e) = c±{e) 


l[Gk{a,e)[Hk{a,e)]^' 
.k=l 


exp 


/ 

.0 



t smh t ZTT 



, (B.l) 



Pfc(2|« 




- 1 




3f 


Pfc(l,^)2Pfc(2 


a 


,0)Pk{2\a 


+ 2,^) 



with C+ = 1, C-{e) = cosh f and 

P, l'9U 



1 + 



(A; + (l + x)/2) 

^/(2vr) 



(2A;+ (2|a| + 3)/2) 
and the exponents are given by 



{k + {l-x)/2) 



fc(fc+i) 

4 



1 + 



(2A;- (2|a| - l)/2) 



:b.2) 



;i - coshi) (1 + cosh((l - 2\a\)t)) Dn{t) ^ cosh((l - 2\a\)t)e 



-Ant 



2sinh^i coshi 

with Dn{t) = [{n + l)(n + 2) - 2n(n + 2) + n{n + l)e-^*]e-2«*. 

If one choose n = 0, this reduces to Eq.(5.26). For the fast convergence, one 
can increase n although the final expression is independent of n. 
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